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^ ■ Abstract 

We consider a Wheeler delayed-choice experiment based on the Mach-Zehnder Interferometer. 
Since the development of the causal interpretation of relativistic boson fields there have not been 
any applications for which the equations of motion for the field have been solved explicitly. Here, 
we provide perhaps the first application of the causal interpretation of boson fields for which the 
equations of motion are solved. Specifically, we consider the electromagnetic field. Solving the 
equations of motion allows us to develop a relativistic causal model of the Wheeler delayed-choice 
, Mach-Zehnder Interferometer. We show explicitly that a photon splits at a beam splitter. We 

also demonstrate the inherent nonlocal nature of a relativistic quantum field. This is particularly 
revealed in a which-path measurement where a quantum is nonlocally absorbed from both arms 
of the interferometer. This feature explains how when a photon is split by a beam splitter it 
nevertheless registers on a detector in one arm of the interferometer. Bohm et al £Q have argued 
that a causal model of a Wheeler delayed-choice experiment avoids the paradox of creating or 
changing history, but they did not provide the details of such a model. The relativistic causal 
model we develop here serves as a detailed example which demonstrates this point, though our 
^J*. model is in terms of a field picture rather than the particle picture of the Bohm-de Broglie 

nonrelativistic causal interpretation. 

^ 1 INTRODUCTION 

In 1978 Wheeler |2] described seven delayed-choice experiments. The experiments are such that the 
choice of which complementary variable to measure is left to the last instant, long after the relevant 
interaction has taken place. Of the seven experiments the delayed-choice experiment based on the 
Mach-Zehnder interferometer is the simplest for detailed mathematical analysis. Here we present 
a detailed model of this experiment based on the causal interpretation of the electromagnetic field 
(hereafter referred to as CIEM) CIEM is a specific case of the causal interpretation of boson 
fields. The experimental arrangement of the delayed-choice Mach-Zehnder interferometer is shown in 
figure 

A single photon 1 enters the interferometer at the first beam splitter BS\. The two beams that 
emerge are recombined at the second beam splitter BS2 by use of the two mirrors Mi and M2. C and 
D are two detectors which can be swung either behind or in front of BS2- The detectors in positions 
C\ and D\ in front of BS2 measure which path the photon traveled and hence a particle description 
is appropriate according to the orthodox interpretation. With the counters in positions C2 and D2 
after BS2 interference is observed and a wave picture is appropriate. A phase shifter P producing a 
phase shift <fi is placed in the /3-beam. For = and a perfectly symmetrical alignment of the beam 
splitters and mirrors, the cZ-beam is extinguished by interference and only the c-beam emerges. 

'email address: pan.kaloyerou@wolfson.ox.ac.uk 

1 Here we use the term 'photon' to mean a single quantum of energy of the electromagnetic field. We do not imply 
any particle-like properties by this term. We have shown in reference 0] that a single photon state is a nonlocal plane 
wave spread over space. 



1 




Fig.l. Delayed-choice Mach-Zehnder interferometer 

If we attribute physical reality to complementary concepts such as wave and particle concepts, then 
we are forced to conclude either that (1) the history of the micro-system leading to the measurement 
is altered by the choice of measurement, or (2) the history of the micro-system is created at the time 
of measurement. 

Wheeler [2] [5], following Heisenberg UJ, in some sense attributed reality to complementary con- 
cepts following measurement and adopted view (2) above, namely, that history is created at the time 
of measurement. Thus he states, 'No phenomenon is a phenomenon until it is an observed phe- 
nomenon,' [7j. He adds that 'Registering equipment operating in the here and now has an undeniable 
part in bringing about that which appears to have happened' p|. Wheeler concludes, 'There is a 
strange sense in which this is a "participatory universe" ' [S]. 

In Wheeler's description, the question of the possibility of creating a causal paradox is raised. 
One can argue, however, in the spirit of Bohr, that Wheeler-delayed choice experiments are mutually 
exclusive in the sense that if the history of a system is fixed by one experiment, this history cannot 
be affected by another Wheeler delayed-choice experiment. But, it is not obvious that the paradox 
can be avoided in this way. 

Bohr and Wheeler share the view that 'no phenomenon is a phenomenon until it is an observed 
phenomenon' but Bohr differs from Wheeler (and Heisenberg) in that he denies the reality of com- 
plementary concepts such as the wave concept and the particle concept. We summerize the features 
of Bohr's principle of complementarity [01513531 as follows: (1) Pairs of complementary concepts 
require mutually exclusive experimental configurations for their definition, (2) Classical concepts are 
essential as abstractions to aid thought and to communicate the results of experiment, but, physical 
reality cannot be attributed to such classical concepts, and (3) The experimental arrangement must 
be viewed as a whole, not further analyzable. Indeed, Bohr defines "phenomenon" to include the 
experimental arrangement. Hence, according to Bohr a description of underlying physical reality is 
impossible. It follows from this that the complementary histories leading to a measurement have no 
more reality than the complementary concepts to which the histories are associated. According to 
Bohr, then, complementary histories, like complementary concepts, are abstractions to aid thought. 

In fact, Bohr had anticipated delayed-choice experiments and writes |12j . '...it obviously can make 
no difference as regards observable effects obtainable by a definite experimental arrangement, whether 
our plans of constructing or handling the instruments are fixed beforehand or whether we prefer to 
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postpone the completion of our planing until a later moment when the particle is already on its way 
from one instrument to another'. Bohr also considers a Mach-Zehnder arrangement but not in 
the delayed-choice configuration. 

Complementarity is not tied to the mathematical formalism. Jammer writes |14| . 'That comple- 
mentarity and Heisenberg-indeterminacy are certainly not synonymous follows from the simple fact 
that the latter... is an immediate mathematical consequence of the formalism of quantum mechanics 
or, more precisely, of the Dirac-Jordan transformation theory, whereas complementarity is an extra- 
neous interpretative addition to it'. Indeed, the whole process from the photon entering BS\ to the 
final act of measurement is described uniquely by the wavefunction (or the wave functional if quantum 
field theory is used, as we shall see). The mathematical description leading up to the measurement 
is completely independent of the last instant choice of what to measure. The wave function or wave 
functional develops causally. Indeed, it is because in the Causal Interpretation mathematical elements 
associated with the wave function or wave functional are interpreted directly that a causal description 
is possible. Wheeler's assertion that a present measurement can affect the past is seen not to be a 
consequence of the quantum formalism, but rather, rests on an extraneous interpretative addition. 
The Bohr view can also be criticized. The denial of the possibility of a description of underlying 
physical reality seems a high price to pay to achieve consistency. 

Clearly, in a causal model of the delayed-choice experiment the issue of changing or creating 
history is avoided. The history leading to measurement is unique and completely independent of the 
last instant choice of what to measure. There is no question of a present measurement affecting the 
past. Bohm et al pQ provided just such a causal description of the Mach-Zehnder Wheeler delayed- 
choice experiment based on the Bohm-de Broglie causal interpretation ^5] though in general 
terms without solving the equations of motion. In this nonrelativistic model, electrons, protons etc. 
are viewed as particles guided by two real fields that codetermine each other. These are the R and 
S'-fields determined by the wave function ip(x, t) = R(x, t) exp[z5(a;, t)/K\. The particle travels along 
one path which is revealed by a which-path measurement (detectors in front of -BS2). The R and 
S'-fields explain interference when the detectors are positioned after BS2- 

Attempts to extend the Bohm-de Broglie causal interpretation to include relativity led to the 
causal interpretation of Boson fields HZ| CD H§| G3 °f wn ich CIEM is a particular example. In CIEM 
the beable is a field; there are no particles. Here, we apply CIEM to the Wheeler delayed-choice 
Mach-Zehnder interferometer. In particular, we set up and solve the equations of motion for the field. 
The solutions allow us to build a detailed causal model of the experiment. In CIEM the basic ontology 
is that of a field, not of a particle as in the Bohm-de Broglie nonrelativistic causal interpretation. 
We will see that a quantum field behaves much like a classical field in many respects but not in 
all. The essential difference is that a quantum field is inherently nonlocal. We will show explicitly 
that a photon is split by a beam splitter. In this case, we will have to show how in a which-path 
measurement, despite being split by the beam splitter, a photon registers in only one of the detectors. 
We will do this by modeling the detectors as hydrogen atoms undergoing the photoelectric effect, and 
show, using standard perturbation theory, that a photon is absorbed nonlocally from both beams by 
only one of the atoms. Since we have a wave model interference is explained in the obvious way. In 
the next section we will briefly summerize CIEM, and in section 13 we apply CIEM to the Wheeler 
delayed-choice Mach-Zehnder interferometer. 



In what follows we use the radiation gauge in which the divergence of the vector potential is zero 
W.A(x,t) = 0, and the scalar potential is zero (/>(x,t) = 0. In this gauge the electromagnetic field 
has only two transverse components. Heavyside-lorentz units are used throughout. 

Second quantization is effected by treating the field A(x, t) and its conjugate momentum II{x, t) 
as operators satisfying the equal-time commutation relations. This procedure is equivalent to intro- 
ducing a field Schrodinger equation 



2 OUTLINE OF CIEM 




(i) 
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where the Hamiltonian density operator H. is obtained from the classical Hamiltonian density of the 
electromagnetic field, 

H = \{E 2 + B 2 ) = \[c 2 n 2 + (V x A) 2 }, (2) 

by the operator replacement 77 — > — iUS/SA. A' is shorthand for A(x',t) and S denotes the 
variational derivative 2 . The solution of the field Schrodingcr equation is the wave functional <&[A, t]. 
The square of the modulus of the wave functional |^[A,i]| 2 gives the probability density for a given 
field configuration A(x, t). This suggests that we take A(x, t) as a beable. Thus, as we have already 
said, the basic ontology is that of a field; there are no photon particles. 

We substitute <P = R[A, t] exp(iS[A, t]/h), where R{A, t] and S[A, t] are two real functional which 
codetermine one another, into the field Schrodinger equation. Then, differentiating, rearranging and 
equating imaginary terms gives a continuity equation: 

dR 2 2 f 5 / 2 5S \ , , n 
+ c / t-tt R -m dx' = 0. 



dt J 5 A' V 6 A 

The continuity equation is interpreted as expressing conservation of probability in function space. 
Equating real terms gives a Hamilton- Jacobi type equation: 

dS 1 ffSS\ 2 o ,„ „ A2 ( h 2 c 2 S 2 R\ , , n 

m + 2j{sA>) c + ( VxA) + {-—jA^) dx = - (3) 

This Hamilton- Jacobi equation differs from its classical counterpart by the extra classical term 

1 f h 2 c 2 S 2 R 



f h'c z 5 Z R , . 
■J-R-5A^ dX > 



which we call the field quantum potential. 

By analogy with classical Hamilton- Jacobi theory we define the total energy and momentum 
conjugate to the field as 

E = _dS[A] n _SS[A] 



dt ' 5A 

In addition to the beables A(x, t) and JJ(x, t) we can define other field beables: the electric field, 
the magnetic induction, the energy and energy density, the momentum and momentum density, the 
intensity, etc. Formulae for these beables are obtained by replacing II by SS/5A in the classical 
formula. 

Thus, we can picture an electromagnetic field as a field in the classical sense, but with the ad- 
ditional property of nonlocality. That the field is inherently nonlocal, meaning that an interaction 
at one point in the field instantaneously influences the field at all other points, can be seen in two 
ways: First, by using Euler's method of finite differences a functional can be approximated as a 
function of infinitely many variables: <&[A, t] — > <P(Aj, A%, . . . , t). Comparison with a many-body 
wavefunction ip(xi,x 2 , ■■■,t) reveals the nonlocality. The second way is from the equation of motion 
of A(x,t), i.e., the free field wave equation. This is obtained by taking the functional derivative of 
the Hamilton- Jacobi equation, 

V 2 4-l^^^ 
c 2 dt 2 SA' 

In general SQ/SA will involve an integral over space in which the integrand contains A(x,t). This 
means that the way that A(x,t) changes with time at one point depends on A(x,t) at all other 
points, hence the inherent nonlocality. 



2 For a scalar function <f> the variational or functional derivative is denned as = A — I f a<t> \ ' ^H- For a 
vector function A we have defined it to be = j^—i + j^— j + j^— k. 
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2.1 Normal mode coordinates 



To proceed it is mathematically easier to expand A(x. t) and II(x, t) as a Fourier series 

A^tJ^^^W^, n(x,t) = -^Y i i klt * klt (t)e- ih -«, (4) 

y 2 — n 



where the field is assumed to be enclosed in a large volume V = L 3 . The wavenumber k runs from 
— oo to +oo and /i = 1, 2 is the polarization index. For A (a;, t) to be a real function we must have 



(5) 



Substituting eq.'s J2J an d @ into eq. Q gives the Schrodingier equation in terms of the normal 
modes q k ^. 



2„2 



d 2 >P 



at' 



(6) 



The solution ^(q k n, t) is an ordinary function of all the normal mode coordinates and this simplifies 
proceedings. We substitute <P = R(q kfl ,t)exp[iS(q klJr ,t)/Ti\, where R(q k ^,t) and S(q kfJi ,t) are real 
functions which codetermine one another, into eq. 10 • Then, differentiating, rearranging and equating 
real terms gives the continuity equation in terms of normal modes: 



OR 2 
dt 



E 



if 



OS 



d 



R 2 



dS 



= 0. 



2 dq k ^ \ dq* k J 2 dq^ \ 8q k ^ 
Equating imaginary terms gives the Hamilton- Jacobi equation in terms of normal modes: 



W + E 



<f__dS_^S_ k 2 , 



d 2 R 



2R dqfdq 



= 0. 



The term 



E 

k ft 



h 2 c 2 d 2 R 

2R dql^q^ 



(7) 



(8) 



is the field quantum potential. Again, by analogy with classical Hamilton- Jacobi theory we define 
the total energy and the conjugate momenta as 



F- 9S 



TTfeu = 



as 

dq k fj, 



dS 



The square of the modulus of the wave function \<P(q kfl) t)\ 2 is the probability density for each q kfl (t) 
to take a particular value at time t. Substituting a particular set of values of q k ^{t) at time t into eq. 
(|4*|) gives a particular field configuration at time t, as before. Substituting the initial values of q k ^(t) 
gives the initial field configuration. 

The normalized ground state solution of the Schrodinger equation is given by 

with TV = Tl'kLiik/ficir)^ 3 . Higher excited states are obtained by the action of the creation operator 



\fnk~t- 



3 The normalization factor N is found by substituting gjf = f k)i + ig kil and its conjugate into #o and using the 
normalization condition \$o\ 2 df kll dg kli = 1, with df kfl = df kll df kl2 df k2l . . ., and similarly for dg k/1 . 
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For a normalized ground state, the higher excited states remain normalized. For ease of writing we 
will not include the normalization factor N in most expressions, but normalization of states will be 
assumed when calculating expectation values. 

Again, the formula for the field beables are obtained by replacing the conjugate momenta Tr k fi 
and 7r£ M by dS/dq kl i and dS/dq k ^ in the corresponding classical formula. The following is a list of 
formulae for the beables: 



The vector potential A(x,t) is given in eq. Q. The electric field is 



E(x,t) = -c II(x,t) 
The magnetic induction is 



1 dA 

c~dt 



c ^ „ dS 



-ih.x 



B(x, t) = V x A(x, t) = — T x ikn)qkn(i)e 



ik.x 



(9) 



(10) 



kfj. 



We may also define the energy density, which includes the quantum potential density (sec reference 
\21\). but we will not write these here as we will not need them. The total energy is found by 
integrating the energy density over V to get, 



F- 9S V 



kfi 

The intensity is equal to momentum density multiplied by c 2 

„2 



Ti 2 c 2 d 2 R 
2R dq* kli dq kll 



i( Xl t) = c 2 g = ^- 



EE 

kfi k' fi' 



i k ^> x (fe x i kfl ) 



dS 
dquu 



-qku,e 



i(k — k ).x 



(11) 



We have adopted the classical definition of intensity in which the intensity is equal to the Poynting 
vector (in heavyside-lorentz units), i.e., J = c(E x B). The definition leads to a moderately simple 
formula for the intensity beable. We note that the definition above contains a zero point intensity. 
But, because J is a vector (whereas energy is not) the contributions to the zero point intensity from 
individual waves with wave vector k cancel each other because of symmetry; for each k there is 
another k pointing in the opposite direction. The above, however, is not the definition normally used 
in quantum optics. This is probably because although it leads to a simple formula for the intensity 
beable it leads to a very cumbersome expression for the intensity operator in terms of the creation 
and annihilation operators: 



-he 2 



EE 

kfi k' fi' 



k k' 

-n£kfi x (k! x e k 'fi') - — (k x i k fi) x ik'fi' 



In quantum optics the intensity operator is defined instead as I = c( x B~ — B~ x 
leads to a much simpler expression in terms of creation and annihilation operators 



(12) 



tic 2 
~V~ 



kfj, k' yJ 



E+), and 



(13) 



This definition is justified because it is proportional to the dominant term in the interaction Hamil- 
tonian for the photoelectric effect upon which instruments to measure intensity are based. We note 
that the two forms of the intensity operator lead to identical expectation values and perhaps further 
justifies the simpler definition of the intensity operator. 

From the above we see that objects such as q k fi, ^kfi, etc. regarded as time independent operators 
in the Schrodinger picture of the usual interpretation become functions of time in CIEM. 
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For a given state <P(qk^,t) of the field we determine the beables by first finding dS/dqk^ and its 
complex conjugate using the formula 



This gives the beables as functions of the qk^{t) and 5jL(i)- The beables can then be obtained in 
terms of the initial values by solving the equations of motion for q%At). There are two alternative 
but equivalent forms of the equations of motion. The first follows from the classical formula 

dC 1 dq* k ^ 



c2 dt ' 

where C is the Lagrangian density of the electromagnetic field, by replacing 717.^ by dS/dqu^,. This 
gives the equation of motion 

1 dq^jt) _ dS 

c 2 dt dq kfl (ty 1 ' 

The second form of the equations of motion for q kll is obtained by differentiating the Hamilton 
Jacobi equation Q by gj£ . This gives the wave equations 

1 dS*^ 2 „ 8Q 

+ K( iku = --a — • ( 16 ) 



C 2 dt 2 ^ dq k ^ 

The corresponding equations for are the complex conjugates of the above. These equations of 
motion differ from the classical free field wave equation by the derivative of the quantum potential. 
From this it follows that where the quantum potential is zero or small the quantum field behaves like 
a classical field. In applications we will obviously choose to solve the simpler eq. lJT5)l. 

In the next section we apply CIEM to the Mach-Zehnder Wheeler delayed-choice experiment. 



3 A CAUSAL MODEL OF THE MACH-ZEHNDER WHEELER 
DELAYED-CHOICE EXPERIMENT 

Consider the Mach-Zehnder arrangement shown in figure 2] BS\ and BS2 are beam splitters, Mi 
and M2 are mirrors and P is a phase shifter that shifts the phase of a wave by an amount 4>. In what 
follows we will assume for simplicity that the beam suffers a ir/2 phase shift at each reflection and a 
zero phase shift upon transmission through a beam splitter. In general, phase shifts upon reflection 
and transmission may be more complicated than this. The requirement is that the commutation 
relations must be preserved. The latter is a more stringent requirement than energy conservation (or, 
equivalently, of photon number conservation) alone |22| . The polarization unit vector is unchanged 
by either reflection or transmission. 



3.1 Input Region 

The input region is the region before the first beam splitter BS\. The incoming beam is a Fock state 
containing one quantum: 



From the formula (|I4|I we can find the S corresponding to the state <Pi (q^ , t) . Using this result in 
eq. Ijl 5(1 gives the equations of motion: 

^ _ cl ^ _ ^ _ ^ _ „ for k ±ko 
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The solutions are easily found: 



< M oW=9oe^ ot+eo) , q*^(t) = q^oe^ for k ± ±k . (18) 

(jo and gfc^o are constant amplitudes, and 9o an d Cfc/xo are constant phases all fixed at t = 0. These 
initial conditions cannot of course be precisely determined but are given with some probability found 
from \$ i {...ql li ...,...q kll ,t)\ 2 . 

We find expressions for the beables in the input region by substituting the dS/dq ko ^ Q and dS/dq kfl 
found from eq. I|17l) together with solutions l|18|l into the formulae for the beables. In finding the 
expression for the beables we use formula (JSJ). Defining 9 = k .x — uj t — 9 the beables are: 

a, \ 2 _ g A (x) 

A{x,t) = — r£fc oAto gocose + , , 

1/2 V 2 

E(x,t) = -^QoEfcoMO sin0 O, (19) 

—2 o fa;) 

B(x,t) = — g (fco x £fc 0M0 )sinO + B ! , 

1/2 V' 2 

with 

3a( :e )= X! i k^e ik - x , g B {x) = V x g A {x) =i ^ (fc x i k ^)q k ^e lk - x . 

kfj, kfj, 



3.2 Region I 

We consider the state <Pj in region I and determine from this state the corresponding beables. Region 
I is the region after the phase shifter P and the mirror Mi but before -BS^- 

The incoming beam <Pj is split at BS\ into two beams: the a and /3-beams 4 . The a-beam 
undergoes a 7r/2 phase shift at M\ and becomes $„e"/ 2 = i$> a . The /3-beam undergoes two n/2 
phase shifts followed by a <fi phase shift and becomes ^e^e 1 * = —'Ppe 1 ^'. Also multiplying by a l/\/2 
normalization factor the state <Pi in region I becomes 

#j = (t«P Q - ^e*) , (20) 

where ^ Q and <£g are solutions of the normal mode Schrodinger equation © and are given by 

*«(«*„.*) = (^) 2 <,^oe-^ ct , ^faW) - 2 /^<2>oe-^ ct . 

The magnitudes of the k- vectors are equal, i.e., k a = kp = fco. 
By using formula l|14|l to first determine 5, 



2i 



ikct - 2ik ct + In (ia% afla - P^^e 1 ^ - In (-ia kafla - (3 ijw e 



(21) 



we can determine dS/da kaila , dS/d/3kptt , and dS/dq^. Substituting the latter into eq. 1|15[1 gives 
the equations of motion for M , (3% ^ and q^: 

da l u 2 9S -tic 2 1 

' - c 2 i^— = —- izjjp (22) 



dt da kafla 2 [ia ka ^ a + Pk^^e' 



4 Some workers insist that two inputs into the beam splitter must be used even when one of the inputs is the 
vacuum 1231 . while other workers use a single input 1241 1251 . In passing, we mention that Caves 1261 . in connection 
with the search for gravitational waves using a Michclson interferometer, suggests, as one of two possible explanations, 
that vacuum fluctuations are responsible for the 'standard quantum limit' which places a limit on the accuracy of 
any measurement of the position of a free mass. We will use the single input approach for BSi as it simplifies the 
mathematical analysis while all essential results remain the same. 
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dt 




2i [ia ka „ a + pk^e-rt] ' 



(23) 



0, for k 7^ ±/c Q , ±kp. 



(24) 



dt 



dqkfi 



Eq.'s (|22J) and (|23|) are coupled differential equations. The coupling shows that the a and /3-beams 
depend nonlocally on each other. Taking the ratio of Eq.'s 11221 and H23|) gives the relation 



which can be used to decouple the two differential equations. The decoupled differential equations 
and eq. (|24|l are easily solved to give 



= aoe l( ^ t+CTo) , Pt^(t) = A)e^ t+To) , g£„(t) = q k ^e<^ for k + ±k a ,±kp, (26) 



where oo and To are integration constants corresponding to the initial phases, and ao and (3q are 
constant amplitudes. Different values of the constants qk^o and Ck^o in different regions of the 
interferometer would correspond to amplitude and phase changes of the normal mode coordinates of 
the ground state. These changes to the normal mode coordinates of the ground state, if indeed they 
occur, do not lead to any observable differences in measured physical quantities so that there is no 
way (at present) to detect them. For this reason we choose qk^o and (kfio to have the same values in 
all the regions of the interferometer. The omega's, uj a = Tic 2 /Aa^ and u>p = hc 2 /A[3 2 , are nonclassical 
frequencies which depend on the amplitudes ao and Po- 

Substitution eq.'s l12(:i[| into eq. (|25|l gives the following relations among the constants associated 
with a fectjUa (t) and f3 k0 ^(t): 



Substituting eq.'s l|27|) into io a or ujp shows that uj a = top. 
3.3 The beables in region I 

In this section we obtain explicite expressions for the beables A(x,t), E(x,t), B(x,t) and I(x,t) as 
functions of time and the initial values. The expression for the energy density is very cumbersome 
and is not as useful in the present context as the intensity. For this reason we will not give the energy 
density here. For the same reason we will also leave out the quantum potential density, though we 
will need the quantum potential. We note that in what follows i ko n = ^k a ^ a = Sk^^p, where £fe oAto 
is the polarization of the incoming wave. 

To find the beables as explicite functions of position and time we substitute the partial derivatives 
of S with respect to the normal mode coordinates found in eq.'s l(22|l . 1(23 jl . and Ij24(l together with the 
solutions for ctk at i a (t), Pkpupit) and qk^t) given in eq.'s f2"fi|l into the formulae for the beables given 
in eq.'s (@J, ©, (| 1 0(1 . and l(f f \ . Eq. 10 is used in deriving the beable expressions. After lengthy 
manipulation and simplification, and defining Q a = k a .x — uj a t — erg and Qp = kp.x — ujpt — To we 
get: 



(25) 



O0 = T + <f> + -. 



(27) 




(28) 



Bi(x,t) 



— r [(ka x £fc QM Jaosine Q + {kp x e fe(3/ ^ )/3 sin Qp] + 




Ii(x,t) 



he 2 f 

— — (k a + kp — k a cos 20 Q — kp cos 2<dp) 




with 



u^x) = ^k^qk^e lk - x , vi(x) =V Xui(x) =i Y ( k x £kn)qky.e ik - x , (29) 

k^±k a ,±k B k^±k a ,±k B 

fj(x) = ihc 2 Y £ Xm« x ( fe x £k^)qk^e lk - x , 9i(x,t) = sin6a, +sin6,a. 

What the above beables show, and the point we want to emphasize, is that the single input photon 
is split by the beam splitter BSi into two beams. Each beam carries half the total momentum, 



, f It hk n hkn 

Gi= \ QjdV = / -4dV = ^ + ^ = hk 
'v Jv 



c 2 2 ' 2 



and energy. The energy distribution given by the Hamilton- Jacobi equation (JZJ) turns out to be very 
cumbersome to calculate, but the total energy is easily found from eq. I|21ll : 

__, dS . , \ - Tick 

E = -m =hch + ^—- 

k 

There is no question of the whole photon choosing a single path through the beam splitter. We note 
that since <!>i is an eigenstate of the energy and momentum operators, the expectation values of these 
operators are equal to the energy and momentum beables. 

We see that the field beable behaves much like a classical electromagnetic field, but there are 
two differences: The first difference is that frequencies oj a — u)p have a nonclassical dependence on 
the amplitude of the field. The second and more significant difference is that each beam depends 
nonlocally on the other beam. This is shown, as we have already mentioned, by the fact that the 
equations of motion of the two beams, eq. (|22|l and eq. Ij23(l . are coupled differential equations. That 
is to say, the behaviour of each beam depends nonlocally on the behaviour of the other beam. 

The nonlocal dependence of each beam on the other can also be seen from the wave equations of 
a k a n a and PkBUB- These can be found by inserting the total quantum potential in region I (found by 
using formula JgJ), 

„ 1 \ ^ , \ ^ tick h 2 c 2 

kfi k 

into the wave equation I jlfijl and differentiating. This gives 



[Hi 



1 d 2 a kafla -h 2 c 2 [a kalla -if3k af t 3 e # ] 



c? dt 2 2(h}h I ) 2 



c 2 dt 2 2{h*hi) 2 

with hi = —ictk a ii a — Pk B ii B exp(— i<f>). In each wave equation the right hand side depends on functions 
from both beams and therefore indicates a nonlocal time dependence of each beam on the other. 

Using eq. (|29|l it is easy to show that the above expressions for the A](x, t), Ei(x, t), and Bj(x, t) 
beables satisfy the usual classical relations E = —(l/c)dAi/dt and Bj = V x Aj. 

There are a number of ways to establish a connection between the initial (constant) amplitudes 
and phases of the input region and region I. One convenient way is to compare the electric field 
beables in the a and /3-beams, eq. (|28|) . with those obtained by beginning with the input electric field 
beable, eq. (|19|l . and inserting the appropriate amplitude and phase changes as it splits at BSi and 
then passes through Mi, and the phase shifter P into region I. This comparison gives the relation 
between the initial constants in the input region and in region I: 

olq = Po = -j=, o- =e -— , t = O - cp - IT. 
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3.4 Which-path measurement 

From the above we see that the beam is split into two parts. With the detectors positioned before BS2 
we know that only one detector will register the absorption of a quantum of electromagnetic energy. 
We must now show how this comes about even though the incoming quantum of electromagnetic 
energy is split into two parts. 

To do this we consider an idealized measurement using the photoelectric effect for a position 
measurement. We model the detectors in positions C\ and D\ in figure ^ as hydrogen atoms in their 
ground state. We assume that the incoming electromagnetic quantum has sufficient energy to ionize 
a hydrogen atom. Each beam interacts with one of the hydrogen atoms. To see what will happen we 
will focus on the interaction of the a-beam with the hydrogen atom at position C\. 

We will treat the hydrogen atom nonrelativistically and picture it as made up of a proton and an 
electron particle according to the Bohm-de Broglie nonrelativistic ontology. From the perspective of 
the description of the electromagnetic field, this nonrelativistic approximation of the atom compared 
to a relativistic treatment will involve only a minor differences in accuracy, but will not alter in anyway 
the model of the field in the interaction process. In the authors opinion, a satisfactory relativistic 
fcrmion ontology has not yet been achieved. A fully relativistic treatment may therefore involve a 
change in the ontology of fermions that we have assumed here. 

The interaction of the electromagnetic field with a hydrogen atom is described by the Schrodinger 
equation 

ih— = (H R + H A + Hj)®. 

Hr, Ha and Hj are the free radiation, free atomic, and interaction Hamiltonians, respectively, and 
are given by 

Hr = J2 («t,«* + I) H a = ^V 2 + V(x), Hj^^^yY: ^e*-^.V, 

with LUk — kc and /i = m e m n / '(m e + m n ) is the reduced mass. The interaction Hamiltonian is derived, 
as usual, using the Pauli minimal coupling. Only the first order terms are retained. We also drop the 
term in Hi containing the creation operator since a hydrogen atom in its ground state cannot emit a 
photon. 

The initial combined state of the radiation and atom is 

<iVi(« fc „,s,t) - ^(g feA1 ) Uj ( : r)e-^ ct e-S fc ^/ 2 e -^*A. 

We take the initial state of the radiation 4>i h to be eq. Ij20|l with the subscript I replaced by 7^. 
Th initial state of the hydrogen atom is taken to be its ground state, 



u l {x 1 t) = -L=e r l a e^ E " t l h 1 



where a = 4nh 2 / /j,e 2 is the Bohr magneton. 
We assume a solution of the form 

<P(q k „X,t) = EE^l.W + fl S,n(*)]^,(ft.)«n(^)e" iEwt/a ) (30) 

where we have retained only the zeroth and first order expansion coefficients, and where En = 
E^ k + E en . Ejf k is the energy of the field including the zero point energy Eq. E en = Tik$c — E e i is 
the kinetic energy of the liberated electron and E e i is the ionization energy of the atom. $N k is the 
final number (Fock) state, and u n (x) — exp(ik en .x) /VV is the outgoing normalized electron plane 
wave. To find the expansion coefficients we use the formulae below derived using standard perturbation 
theory: 

a ivL«(*) = S N hf >i h Jni, (31) 
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,(!) 



iiJM, n I, it ji 
2 "hn n >'hii' dt, 



tiUNk»n,i hft .i = E Nk ^ n - E Ihlii = E Nkti + E en - E Ikfi - E ei = E Nk ^ n j k ^. 
Evaluating the time integral in eq. 1)32(1 gives 



1 

if) 



g "k^-'kn' (It = 



E 



Nk^.n,I kll i 



(32) 

(33) 
(34) 



(35) 



The modulus squared of eq. I|35|l in the limit t — > oo becomes 5(EN k ^ n j kfi i)2Trt/h. This corresponds 
to energy conservation. We require the time t over which the integral is taken to be very much longer 
than h/E n , but sufficiently short that a^(t) does not change very much. In this case, it is a good 
approximation to take t — > oo as we have done above. 

After substituting for Hi in eq. I|33|) the matrix element becomes 



H 



N k)1 n,I kli i 



—e / he 



X! ~jf / ®*NkM^)ak^i k Mkn) d 1kti 

kfl V 



The first integral is 



u* n {x)e ik - x i kll .{-iTiV)u l {x) dx 

(36) 



■^=(i- e^SN^oSkkJfifio- (37) 



Using the dipole approximation e lh x = 1 and the fact that the delta functions Skk 5nno in the first 
integral produce ffc M — ► £k n , the second integral of eq. i|36|) is evaluated to give 



Un(x)i koflo .(-iKV)ui(x) dx = 



f^kofj-a -ken 



(l+« 2 fce 2 J 2 



(38) 



where £k ^o -k en is the component of the outgoing electron wave vector in the direction of polarization 
iko/io °f the incident electromagnetic field. 

We draw attention to the fact that the integral (|37|l demonstrates that if interaction takes place at 
all then the atom must absorb one entire quantum and leave the field in its ground state. Substituting 
eq.'s (fTTTT) . (J221, <|3"5f . (jSEJ, EH and ll3"%|l into the assumed solution eq. (|3T))l gives the final solution 



^ ^ i i\ ^o(gV>*)v^ us- i 1 

$ = $I k Alk^ X , t) Y 2^ mn{t)£k aM - k en-^t 



t(fe e 



with 



VOn(t) 




(i - e i 



8?ra 3 



(l + a 2 fc2 n )- 



c-B e „t/B.) 



o lJ,t - 1 k u ' 



Eo nt i kfl i 



(39) 



Eon,l kfl i is given by eq. (£2} with N kll = 0. 

In the solution (|39|l the first term is the initial state and corresponds to no interaction taking 
place. Recall that for a single atom the probability of interaction with the electromagnetic field is 
small. The second term shows that if the interaction takes place then one entire quantum must be 
absorbed. This means that the field energy from both beams is absorbed by only one of the hydrogen 
atoms. Since the arms of the interferometer can be of arbitrary length, and since the duration of the 
interaction with the atom is small, the absorption of the electromagnetic quantum occurs nonlocally. 
This is the second way in which nonlocality enters into the description, and further emphasizes the 
difference from a classical field. In this way, we show how despite the fact that a single quantum is 
necessarily split by the first beam splitter, only one of the counters placed before BS2 fires. 
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3.5 Region II 

With the counters in the Ci and D2 positions interference is observed. Since we are modeling the 
quantum mechanical electromagnetic field as a real field, interference is explained in the obvious way 
even though the field has additional nonclassical properties. Here we want to show the interference 
of the beables explicitely by solving the equations of motion for the normal mode coordinates. 

At BS2 the a and /3-beams are split forming the c and d-beams by interference (see figure We 
call the region after BS2 region II. We want to find the state $u of the field in this region and to 
determine the corresponding beables. 

The part of the a-beam reflected at BS2 suffers another ir/2 phase shift. The part of the /3-beam 
transmitted at BS% suffers no phase shift. Each beam is again multiplied by 1/V2 because of the 
50% intensity reduction at BS2- The two beams interfere to form the c-beam represented by the state 
-(l/2)* c (l + e**). 

The transmitted part of the a-beam experiences no phase change, while the part of the (3 beam 
that is reflected at BS2 suffers a further 7r/2 phase change. These two beams interfere to form the 
d-beam represented by the state (i/2)<P d (l — e*^). Adding these two states gives the state of the field 
in region II 

*iJ = -J*c(l + e i *) + i# ( j(l-e i *). 



2' 



<I C and <S>a are Fock states given by 



Note that the magnitudes of the A:- vectors are unchanged by interaction with optical elements, i.e., 
k c = kd = k a — kp = k . As before, by using formula i|14|) to determine S, 

ikct - 2ik ct + In [-(1 + e^)c^ c + i(l - e**)^ J 

k 

-(l + e-^K^-ia-e" 4 *)^]}, 

we can determine the derivatives of S with respect to the normal mode coordinates. Substituting the 
latter into eq. I|15f) gives the equations of motion of the normal modes: 



dc* kalic 


2 ds 


c 2 h (1 + e~ 


i<j>\ 




dt 


dc kafic 


2 [-t(l + e-^)c kcllc + 


(1- 




d d*k u 

kdf-d 


2 dS 


c 2 h i(l — e 






dt 


ddk df i d 


2 H(l + e-^)c fccMc + 


(1- 


- e~^)d kdt2d 


dt 


2 dS 
dq kfl 


0, for k 7^ ±fe c , ±kd 







(40) 

By taking the ratio of the two coupled differential equations to get the relation 

i(l-e-^)4 oMo = (l + e-^K dW , (41) 

and by following similar steps to those of section 1|3.2(1 , we can find the solutions of the equations of 
motion (gUJ: 

c)UW = c e^ t+xo) , d* killi {t) = d Q e^ t+ Z°\ q* kfl (t) = ft^e**" for k ± ±k c ,±k d , (42) 

with uj c = [7ic 2 (l + cos0)]/4cq and u>d — [?ic 2 (l — cos0]/4iig. The Cq, do and q k u,o are constant 
amplitudes, and xcb £0 an d the £ku,o are constant phases. 

Setting t = in eq. 's l|41|) , and (|42|l and solving the resulting equations gives the relation between 
the initial (constant) values in region II: 

d = -c tan -, xo = fo- (43) 

Substituting eq.'s 1)43(1 into ui c or u>d shows that ui c = uj d . 
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3.6 The beables in region II 

As in section 13.81 to get expressions for the beables we substitute the derivatives of S with respect 
to the normal mode coordinates given in eq.'s (|4C)|1 together with the solutions for the normal mode 
coordinates, eq.'s ((42(1 . into the formulae for the beables given in subsection 12.11 We note that the 
polarization remains unchanged in interactions with optical elements, i.e., £fc cMc = £k d fj, d — £fc Q ^ Q — 
£kpiip — £k no and that we have used eq. JSJ). After lengthy manipulation and simplification, and 
defining O c = k c .x — uj c t — xo and ©d = kd-x — Lu^t — £o we get: 

. , . 2 . uu{x) 
A n {x,t) = — - (£ fccAtc c cosB c + e kdPLd d cosB d ) H - t — , 

Y 2 V 2 

E H (x,t) = ^ ( ^=^(l + cos^)sine c + ^i(l-cos^)sined ) , 
2V* V c o «o / 

— 2 ^//('^) 
B H (x,t) = — [(fe c x £fc cM Jcosine c + (fc d x e fed/id )d sin 6 d ] H j— , 

V 2 V/ 2 

he 2 

Iu(x, t) — [k c (l + cos 4>) + kd(l — cos <ft) — k c (l + cos 4>) cos 29 c 4- fed(l — cos 4>) cos 20<j)] 



2V 

fn( x )9ii(x,t) 
V 



(44) 



with 



u H (x) = £k fJ ,qk f ,e lkx , v u (x)=i ^ (k x i kf _ l )q kf _ l e lk - x = V x u H {x), (45) 

= £ feoMo x (fe x ik^)qk^e lk - x , g H {x,t) = (1 + cos0)sin0 c + (1 - cos sin 6d. 

kyi±k c ,±k d 

Using cq. (|45|) we can again show that the above expressions for the Ajj(x,t), Ejj(x,t), and 
B u (x,t) beables satisfy the usual classical relations En = —{\/c)dA/dt and B n = Vx A H . 

By tracing the input electric field beable, eq. i(19|) . through all of the optical elements and then 
comparing with the electric field beable obtained from the equations of motion (as we did in section 
13.3(1 we obtain relations among the initial (constant) values in all three regions: 

f— (j) /—(f) <j> 7T 

c = V2a cos-, d = -V2a sin-, Xo = So = - g _ o ' 

c = q a cos -, d = -go sin ^, X o = Co = #o - ^ - 7r. (46) 

Substituting eq. ((46(1 into w c or u>d shows that wo = u c = = — 

The beables above clearly exhibit interference. For example, with a phase shift of (f> = the 
d-beam is extinguished, while for <f> = ir the c-beam is extinguished. For Iu(x,t) and the En(x,t) 
beables this is obvious, but for Au(x,t) and Bjj(x,t) we need to use relation {13} . 

In the experiment what is actually observed is the expectation value of the intensity operator, 
eq. O (or eq.® which gives the same results). This expectation value is equal to the long time 
average (which is equal to the cycle average for a periodic function) of the intensity beable, eq.lJSJl, 
i.e., 

1 [ T " /2 he 2 

^I^cycle = f J T /2 IU dt = Qy^ 1 +C ° S ^ +M 1 "COS </>)], 

with T c = 2ir/u! c . The total momentum and energy of each beam, 

Tick 



G u = 



/ -j" dV = — [fe c (l + cos0) + fc ( ;(l-cos0)] =hko, E=Ucko+y^ 

Jy C IV A 



are again equal to the corresponding expectation values, as we would expect. 
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As with the a and /3-beams the two coupled differential equations in (|4l)|l show that the two beams 
are nonlocally connected. This is further demonstrated by the wave equations 



1 d 2 c kalla -2h 2 c 2 [(1 + cos(j>)c kcflc + s'mcj) d kd 



c 2 dt 2 {h* n h u ) 2 



"IV 

; 2 2 I 



1 d 2 d kdt _ ld _ -2h c 2 [sin c fccAtc + (1 - cos<ft)rf fcdAld ] 



c 2 dt 2 (h* n h n ) 2 
obtained by substituting the quantum potential in region II, 



ill 



1 \ - 2 \ - Aifcc 



hkc h 2 c 2 



k/i k 11 11 



with hjj = — [1 + cxp(— «</))]c/ Cc/ic — i[l — exp(— i<fi)]d kdfJid , into the wave equation (|1 ti|> . 



4 CONCLUSION 

We have provided perhaps one of the first applications of the causal interpretation of relativistic 
boson fields where the equations of motion for the field are solved explicitly. In so doing we have 
been able to provide a detailed relativistic causal model of the Wheeler delayed-choice Mach-Zehnder 
Interferometer, i.e., a description of the physical reality that underlies the experiment. We have shown 
explicitly that a single photon is split by a beam splitter. We have shown that the beables representing 
quantities such as the electric field and the magnetic induction behave much like their classical 
counterparts. One similarity is that the expressions for the electric field and magnetic induction 
beables in terms of the vector potential beable are the same as for their classical counterparts. They 
differ from their classical counterparts in that the beables oscillate with a nonclassical frequency 
which depends on the amplitude of the wave. A more significant difference is the inherent nonlocality 
of a relativistic quantum field. This nonlocality is revealed in two ways: First, because the time 
dependence of a beam in one arm of the interferometer depends nonlocally on the time dependence 
of the beam in the other arm. Second, in a which-path measurement an entire quantum is absorbed 
nonlocally from both arms of the interferometer by a detector placed in one arm. This feature explains 
how when a photon is split by a beam splitter it nevertheless registers in a detector placed in one 
arm of the interferometer. Wheeler concludes from his hypothetical delayed-choice experiments that 
history is created at the time of measurement. Bohr tells us that complementary concepts such 
as wave and particle concepts are abstractions to aid thought, and hence, so also is any historical 
evolution leading to the final experimental result. We have argued that neither Wheeler's nor Bohr's 
conclusions follow from the mathematical formalism of the quantum theory. The Bohm-de Broglie 
nonrelativistic causal interpretation and its relativistic generalization to boson fields demonstrates 
that we neither need to follow Bohr and deny a description of underlying physical reality, nor do we 
need to follow Wheeler and conclude that the present can affect the past. 
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